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GRAVITATIONAL WAVES FROM SUPERMASSIVE BLACK HOLE COALESCENCE IN A HIERARCHICAL 

GALAXY FORMATION MODEL 
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ABSTRACT 

We investigate the expected gravitational wave emission from coalescing supermassive black hole 
(SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs 
in the host galaxies sink to the center of the new merged galaxy and form a binary system. We 
employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering 
scenario to estimate the amplitude of the expected stochastic gravitational wave background owing 
to inspiraling SMBH binaries and bursts owing to the SMBH binary coalescence events. We find that 
the characteristic strain amplitude of the background radiation is h c (f) ~ l(V 16 (//l^Hz)~ 2 / 3 for 
/ < 1/iHz just below the detection limit from measurements of the pulsar timing provided that SMBHs 
coalesce simultaneously when host galaxies merge. The main contribution to the total strain amplitude 
of the background radiation comes from SMBH coalescence events at < z < 1. We also find that a 
future space-based gravitational wave interferometer such as the planned Laser Interferometer Space 
Antenna (LISA) might detect intense gravitational wave bursts associated with coalescence of SMBH 
binaries with total mass M tot < 10 7 M Q at z > 2 at a rate ~ 1.0 yr~ 4 . Our model predicts that burst 
signals with a larger amplitude /iburet ~ 1CV 15 correspond to coalescence events of massive SMBH 
, binary with total mass M to t ~ 10 8 M Q at low redshift z < 1 at a rate ~ 0.1 yr" 1 whereas those with a 

smaller amplitude /iburst ~ 10 -17 correspond to coalescence events of less massive SMBH binary with 
Q\ ' total mass M to t ~ 10 6 Af© at high redshift z > 3. 

Subject headings: black hole physics - galaxies: evolution - galaxies: formation - gravitational waves 
-quasars : general 

o 
o 



o 
o 

(N 



in 



Oh 



1. INTRODUCTION 



In recent years, there has been increasing observational evidence that many nearby galaxies have central supermassive 
black holes (SMBHs) in the mass range of 10 6 — 10 9 M Q , and that their physical properties correlate with those 
of spheroids 5 of their host galaxies. First, the estimated mass of a SMBH in a galactic center, Mbh, is roughly 
O ' proportional to the mass of the s pheroid, M sn h Br nid. The ratio of Mrh /-M^wnM i s 0.001 — 0.006 in each galaxy 
(e.g. iKormendv fc Richstonelll995l iMaeorrian et al.lll998t iMerritt fc Ferrarese 2001b|). Second, the mass of a SMBH 
c/3 ■ correlates with the velocity dispersion of stars in the spheroid, <7 sp heroid, a s Mbh oc cr" phe roid ,n = 3.7 — 5.3 (e.g. 
! iFerrarese fc Merrittl l2000t iGebhardt et~al] 1200(1 IMerritt fc Ferrarese! l207)Tal iTremaine etaLl f2"002L These relations 
suggest that the formation of SMBHs physically links to the formation of spheroids that harbor the SMBHs. 

In order to study the formation and evolution of SMBHs, it is necessary to construct a model that includes galaxy 
formation and merger processes. It has become widely accepted that quasars are fueled by accretion of gas onto 
■ SMBHs in the nuclei of host galaxies. Hydrodynamical simulations have shown that a merger of galaxies drives 
gas fall rapidly onto the center of the merged system and fu els a nuclear starburst leading to spheroid formation 
(e.g. iMihos fc Hernouistlll994 IT996T: iBarnes & Hernauistl H996). Some obser vations of quasar hosts show that many 
quasars reside in spheroids of interacting s ystems or elliptical galaxies (e.g. iBah call ct aLl ll997t iMcLure et al.lll999i 
iDunlop et aTll2003t iKauffmann et alJl2003|) . Thus the major merger of galaxies is a possible mechanism for quasar 
activity, SMBH evolution and spheroid formation. 

In the standard hierarchical structure formation scenario in a cold dark matter (CDM) universe, dark-matter halos 
(dark halos) cluster gravitationally and merge together. In each of merged dark halos, a galaxy is formed as a result of 
radiative gas cooling, star formation, and supernova feedback. Several galaxies in a common dark halo sometimes merge 
together and a more massive galaxy is assembled. When galaxies merge, SMBHs in the centers of the galaxies sink 
toward the center of the new merged galaxy and form a SMBH binary subsequently. If the binary loses enough energy 
and angular momentum, it will evolve to the gravitational wave emitting regime and begin inspiraling, eventually 
coalesces with a gravitational wave burst. 
An ensemble of gravitational waves from a number of inspiraling SMBH binaries at different redshift can be observed 
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as a stochastic b ackground at frequencies ~ In — 1/iHz, which can be detected by pulsar timing measurements 
l)Detweiler II1979J) . Gravit ational wave bursts from c oalescence events can be detected by a Doppler tracking test of 
interplanetary spacecraft l)Thorne fc Braginskvlll97^) . Future space interferometers such as the Laser Interferometer 
Space Antenna LISA might detect quasi-monochromatic wave sources associated with inspiraling SMBH binaries and 
gravitational wave bursts associated with SMBH binary coalescence. 

In order to estimate the feasibility of detecting such gravitational waves, we need to know the SMBH coalescing rate 
that depends on the complex physical processes before coalescence. A possible scenario of forming SMBHs consists of 
four stages: First, a seed massive black hole in the center of a galaxy grows to an intermediate mass black hole (IMBH) 
via run-away collapse in a dense star cluster ijPortegies Zwart. et alJl2004|) . Second, a galaxy that contains an IMBH 
in their center merges with another galaxy with an IMBH. Third, the two IMB Hs sink to the center of gra vitational 
potential well owing to dynamical friction until a "hardening" regime ijVolonteri. Haardt fc Madaul 12003). Finally, 
physical processes such as three-body processes, gas dynamics, the "Koza i " mechanism, etc. drive the binary to the 
gravitational wave emission regime re.g.lBegelman. Blandford fc R.eesll980HMakinoil997tlGould fc RbJ200fllYJ2002 
[ Armitaee fc Nataraianll2002t iBlaes. Lee, fc Socratesll2002t iVolonteri. Haardt fc Madaull2003t iMilosavlievic fc Merrittl 
2003). However , the efficiency to bring two SMBHs together is still unknown. 

To date, a number of attempts have been made to calculate the SMBH coalescence rate taking the physical pro- 
cesses in the "second" and "third" stages into account. S ome authors use phenomenological models of galaxy mergers 
based on number counts of quasars and spheroids fe.g. iTh orne fc Braeinskvl 119761: iFukushige. Ebisuzaki fc Makinol 
1992: H aehneltl 119941: iRai agonal fc Romanil 119911 I.Iaffe fc Backer! 12003) . Others use merger rates of dark ha- 
los (not galaxies) based on forma lism of Press fc Schechterl (119741) and ILacev fc Col3 lfl993T) (e.g. IHaehuelH ITflfll 
iMenou. Haiman fc Narayan an 2001). IVolonteri. Haardt fc Madaul ()2003) calculated the SMBH coalescence rate tak- 
ing conditions under which "sub-halos" w ith SMBHs in a dark halo sink to the dark halo center into account (see also 
IWvithe fc Lo ^ l2003l lSe" sa,na et a,U l200l. However, none of these models include baryonic gas evolution and galaxy 
formation processes. Because SMBH formation process is relevant to spheroids of host galaxies rather than to dark 
halos, we need to evaluate how the baryonic gas processes such as star formation and radiative cooling affect the 
SMBH coalescence rate. 

In this paper, we estimate the SMBH coal escence rate using a new semi-analytic (SA) model of En oki et alJ l)2003l) (an 
extended model of Nagashi ma et alJ l)2001|) ) in which the SMBH formation is incorporated into the galaxy formation. 
Then, we calculate t he spectrum of gravit ational wave background from inspiraling SMBH binaries, based on the 
formulation given bv LTaffe fc B^ckel Ipvil and we compare our result with that from a pulsar timing measurement. 
We also estimate the event rate of gravitational wave bursts from SM BH coalescence events that might be detected 
by future planned space laser interferometers, based on an argument in iTh orne fc Braginskvl l)1976j) . 

In SA models, merging histories of dark halos are realized using a Monte-Carlo algorithm and evolution of baryonic 
components within dark halos is calculated using simple analytic models for gas cooling, star formation, supernova 
feedback, galaxy merging and other processes. SA models have successfully reprod uced a variety of observed features 
of ga l axies, such as their luminosi ty functions, color distri bution s, and so on (e. g. |K^iu^hTiann_et al.lll993t iCole et all 
119941 1200(1 ISomerville fc Primacklll999t iNagashima et alJl200ll l2002j) . iKauffmann fc HaehneltJ (|2000j) introduced a 
unif ied model that includes the formation of both galaxies and SMBHs within the fra mework of their SA model (see 
also lCattane73l2001tlMenci et alJl2003t llslam. Tavlor fc SiHdl2003UGranato et al.l f2004L Our SA model reproduces not 
only these observational featur es but also the ob served present black hole mass function and the quasar luminosity 
functions at different redshifts l)Enoki et al .1 12003 L 

The paper is organized as follows: in section [2] we briefly review our SA model for galaxy formation and SMBH 
growth; in section we calculate the spectrum of gravitational wave background and the event rate of gravitational 
wave bursts; in section 0] we provide summary and conclusions. 

2. GALAXY MERGER / BLACK HOLE COALESCENCE RATE 

In this section we briefly describe our SA model for galaxy formation and the SMBH growth. The details are shown 
in INagashima et alJ (|2001|) and lEnoki et~aT1 1)2003|) . 

2.1. Model of Galaxy Formation 

First, we construct Mo nte Carlo realization s of m erging histories of dark halos from the present to higher red- 
shifts using a method o f ISomerville fc Kolattj ltl"999l). which is based on the extended Press-Schechter formalism 
(jPress fc SchechterlH974l iBond et al.lll99lUBowerlll99ll ILacev fc Coldll993|) . Merging histories of dark halos depend 
on the cosmological model. The adopted cosmological model is a low-density, spatially flat cold dark matter (ACDM) 
universe with the present density parameter, f2 m = 0.3, the cosmological constant, tt\ — 0.7, the Hubble constant 
h = 0.7 (h = i/o/100 km s^ 1 Mpc" 1 ) and the present rms density fluctuation in spheres of 8/i _1 Mpc radius, erg = 0.9. 
Dark halos with circular velocity, V^rc <40 km s , are treated as diffuse accr etion matter. This conditi on comes from 
the estimation of Jeans mass in the ultraviolet background radiation field fe.g. lThoul fc Weinbergl H996'). The adopted 
timestep of merging histories of dark halos is a redshift interval of Az = 0.06(1 + z), corresponding to dynamical 
time scale of dark halos which collapse at redshift z. The highest redshift in each merging path which depends on the 
present dark halo mass, is about z ~ 20 — 30. 

Next, in each merging path of dark halos, we calculate the evolution of the baryonic component from higher redshifts 
to the present. If a dark halo has no progenitor halos, the mass fraction of the gas in the halo is given by f2b/^ m 
where ilb is the baryonic density parameter. Here we use fib = 0.02/i -2 . When a dark halo collapses, the gas in the 
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halo is shock-heated to the virial temperature of the halo (the hot gas). At the same time, the gas in dense regions 
of the halo cools owing to efficient radiative cooling and sinks to the center of the halo and settle into a rotationally 
supported disk until the subsequent collapse of the dark halo. We call this cooled gas the cold gas. Stars are formed 
from the cold gas at a rate of = Af co id /t* , where M co id is the mass of the cold gas and t* is the time scale of 
star formation. We assume r* = T°(Vcirc/300km s -1 )"*. The free parameters of t° and a* are chosen to match the 
observed mass fraction of cold gas in the disks of spiral galaxies. With star formation, supernovae occur and heat up the 
surrounding cold gas to the hot gas phase (supernova feedback). The reheating rate is given by M re heat = /3(Vcirc)M*, 
where /3(T4i rc ) = (Vhot/Vcirc) QhDt - The free parameters of Vhot and ahot are determined by matching the observed local 
luminosity function of galaxies. Given the star formation rate as a function of time, we calculate the luminosity and 
color of each galaxy from the star fo rmation history of the galax y by using a stellar population synthesis model. We 
use the population synthesis code bv lKodama fc Arimoto I l|1997|) . 

When several progenitor halos merge, a newly formed larger dark halo contains at least two or more galaxies which 
originally resided in the individual progenitor halos. We identify the central galaxy in the new common halo with 
the central galaxy contained in the most massive progenitor halos. Note that cooled hot gas accretes to only central 
galaxies. Other galaxies are regarded as satellite galaxies. These galaxies merge by either dynamical friction or random 
collision. Satellite galaxies merge with the central galaxy owing to dynamical friction in the following time scale, 

260 f RaV f V ciic \f M sat \~ x 

where i?H and V c i rc are the radius and the circular velocity of the new common halo, respectively, In A r is t he Coulomb 
logarithm, and M sat is the mass of the satellite galaxy including its dark halo ijBinnev fc Tremaindll987f) . When the 
time passing after a galaxy became a satellite exceeds Tf r i c , the satellite galaxy infalls onto the central galaxy. On the 
other hand, satellite galaxies merge with each other in timescale of random collision. Under the condition that th e 
satellite galaxies gravitationally bound and merge during encounters, the collision time scale is l)Makino fc Hird ll997). 

500 fR H \ 3 f r gal 

Tcoll ~ 



iV 2 \MpcJ V 0.12 Mpc 



-4 / \ 3 

OVal \ / Chalo 



100kms _i 7 \ 300 kms 



T Gyr, (2) 



where N is the number of satellite galaxies, r ga i is a radius of a satellite, and <7 n aio and (7g a i are the ID velocity 
dispersions of the common halo and the satellite galaxy, respectively. With a probability of A£/t co ii, where At is the 
timestep corresponding to the redshift interval Az, a satellite galaxy merges with another randomly picked satellite. 
Let us consider the case that two galaxies of masses mi and ni2(> m±) merge together. If the mass ratio, / = mi/m2, 
is larger than a certain critical value of /bulge (major merger), we assume that a starburst occurs: all of the cold gas 
turns into stars and hot gas which fills the dark halo, and all of the stars populate the bulge of a new galaxy. On the 
other hand, if / < /bulge (minor merger), no starburst occurs and a smaller galaxy is simply absorbed into the disk of a 
larger galaxy. These processes are repeated until the output redshift. We classify galaxies into different morphological 
types according to the -B-band bulge-to-disk luminosity ratio, B/D. In this paper, galaxies with B/D > 2/3, and 
B/D < 2/3 are classified as ellipticals/SOs and spirals, respectively. The parameter /bulge is fixed by a comparison 
with the observed type mix. 

Model parameters are determined by a comparison with observations of galaxies in the local Universe, such as 
luminosity functions and the cold gas mass fraction in spiral galaxies. Our SA model can reproduce galaxy number 
counts and photometric redshift distribution of galaxies in the Hubble Deep Field. The adopted parameters of this 
model are tabulated in table ^ Some of the model pa rameters (as , and fraction of invisible stellar mass T) are 
updated and slightly different from our previous paper ijNagashima et al.ll2001t lEnoki et alJ l2003'). The update of <r 8 
and fib causes only a slight change of T and we have confirmed that the modification hardly affects our results. Using 
this SA model, we can estimate the galaxy merger rate at each redshift. 

2.2. The growth of SMBH 

Let us briefly summarize the growth model of SMBHs introduced bv lEnoki et al.l l)2003j) . In this model, it is assumed 
that SMBHs grow by coalescence when their host galaxies merge and are fueled by accreted cold gas during major 
mergers of galaxies. When the host galaxies merge, pre-existing SMBHs sink to the center of the new merged galaxy 
due to dynamical friction (or other mechanisms such as gas dynamics), evolve to the gravitational wave emission 
regime and eventually coalesce. Although the timescale for this process is unknown, for the sake of simplicity we 
assume that SMBHs instantaneously evolve to the gravitational wave emission regime and coalesce. Gas-dynamical 
simulations have demo nstrated that the major merger of galaxies can drive subs tantial gaseous inflows and trigger 
starburst activity (e.g. Mihos fc Hernauistj fl994. 19961 iBarnes fc Hernauistlll996|) . Thus, it is reasonable to assume 
that during a major merger of galaxies, a certain fraction of the cold gas that is proportional to the total mass of stars 
newly formed at starburst accretes onto the newly formed SMBH and this accretion process leads to a quasar activity. 
Under this assumption, the mass of cold gas accreted on a SMBH is given by 

M acc = / BH AM* iburst) (3) 
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where /bh is a constant and AM^bm-st is the total mass of stars formed at starburst. The free parameter of /bh is 
fixed by matching the obser ved relation between a spheroid mass and a black hole mass MBH/M sp heroid = 0.001 — 0.006 
(e.g. iMaeorrian et al.|l!998|k we find that the favorable value of /bh is nearly 0.03. 

Figure^shows the model prediction of the growth rate of the averaged SMBH mass, (Mbh}- We define the averaged 
SMBH mass growth rate as follows, 

/»> v f Mbh<Pbh(M bh , z)dM BH 

(M B h) = 7-T 7T7 rTTT ' W 

J <Pbh(M B h, z)dM BH 

where Mbh is the mass increase rate of SMBH with A/bh at z, and </>bh(M B h, z) is the black hole mass function. Mbh 
is given by Mbh = (M acc + M coa i)/Ai where M coa i is the mass increment due to SMBH binaries coalescence (the mass 
of smaller SMBH) . This figure shows the SMBH mass growth is mostly due to gas accretion and does not depend on 
initial seed masses. At lower redshift (z < 1), since the cold gas is depleted by star formation, the SMBH coalescence 
becomes a dominant process of the SMBH mass growth. 

Figure|21(a) shows the black hole mass functions in our model at a series of various redshifts. This indicates that the 
number density of the most massive black holes increases monotonically with time in our scenario in which SMBHs 
grow by an accretion of co l d gas and by coalescence. In this figure, we superpose the black hole mass function at z — 
obtained bv iSalucci et alJ l)1999fl . The present mass function in our model is quite consistent with the observation. 
Our galaxy formation model includes dynamical friction and random collision as galaxy merging mechanisms. For 
comparison, in figure [3(b), we also plot the black hole mass functions at z = of other two models: no random 
collision model (no rc model) and no dynamical friction model (no df model). In the no rc model and the no df 
model, mergers owing to random collision and dynamical friction are switched off, respectively. This figure shows 
that the mass function for low mass black holes are determined by random collisions between satellite galaxies and 
that for high mass black holes are influenced by dynamical friction. The shape of black hole mass function depends 
on detailed gas processes. The important contribution of mass increment of SMBHs in central galaxies is the cold 
gas, which is accreted to only central galaxies. SNe feedback removes this cold gas more efficiently in smaller galaxies 
with 14i rc < Vhot = 280km s" 1 (M ga i < 10 12 M Q ). Thus, the growth of the SMBHs in small central galaxies suffers 
from SNe feedback. In the no rc model, SMBHs mainly exist central galaxies. Therefore, the shape of the black hole 
mass function in the no rc model has a bump at high mass end (Mbh ~ 10 9 M Q ) which corresponds to SMBHs in the 
central galaxies (M ga i ~ 10 12 M Q ). On the other hand, in the no df model, high mass SMBHs cannot be produced 
si nce galaxies canno t merge with the massive central galaxy. 

lEnoki et al.l l)2003|) investigated environments and clustering properties of quasars using this SA model which can also 
reprod uce quasar luminosity functions at different redshifts. While our approach is similar to Kauffma nn fe HaehneU 
(2000), their star formation and feedback models are different from ours and their model does not include the random 
collision process. Therefore, their resultant model description differs from ours in equation ©• Using the SA model 
incorporated with this SMBH growth model, we estimate the comoving number density, n c {M\, M2, z)dMidM2<iz, of 
the coalescing SMBH binaries with mass Mi ~ M% + dM\ and M2 ~ M2 + dM-2 at 2 in a given redshift interval dz. 

It is difficult to know how SMBH binaries manage to shrink to a gravitational wave emission regime after their 
host galaxies merge, because all physical processes and conditions related to this problem are still not clear (dynam- 
ical friction, stellar distribution, tripl et SMBH interaction, gas dyn amical effects and so on). Neve rthel ess, many 
authors have tackled this problem fe.g.lBegelman. Blandford fc Reeslll98rtlMakinofll997t iGould fc Rb3l200ft lYull2002t 
lArmitage fc Nataraianll2002tlBlaes. Lee, fc Socratesll2002UMilosavlievic fc Merrittll2003|) . In what follows, we assume 
that SMBHs coalesce simultaneously when host gal axies merge for simplicity. In other words, the efficiency of SMBH 
coalescence is assumed to be maximum. Note that iVolonteri. Haardt fc Madaul l)2003l) constructed a SMBH growth 
model in which the merging history of dark halos with SMBHs in their center is derived from cosmological Monte-Carlo 
simulations. Although they incorporated dynamical evolution of SMBH binaries and triplet SMBH interactions into 
their model, no galaxy formation process is taken into account. 

3. GRAVITATIONAL RADIATION 
3.1. Background radiation from SMBH binaries 

In order to calculate the spectrum of gravitational wave background radiation from SMBH binaries, it is necessary 
to follow the orbit evol ution of each binary. H owever, i t is difficult to set initial c onditions. To avoid this problem, we 
adopt a formulation bv l.Taffe fc Backer! l)2003[) (see also lRaiagopal fc Romanilll995th 

We start calculation of the expected amplitude of gravitational radiation emit ted by a binary. We assume that a 
binary orbit is circular for simplicity 6 . The amplitude is given by llThornelll987l) 

h s (z, /, M U M 2 ) = 4yf (G ^ 5/3 (2nf p f 3 



= 3.5 x 10" 



( M chirp x 


^5/3 


r d(z) - 


-1 


|7(i + *)] 


liO 8 M 0/ 




1 Gpc 




lO" 7 Hz 



2/3 

(5) 



6 Note that the initial orbit of a binary is not necessarily circular. S ince the amplitude and timescale of gravitational wave emissions 
are sensitive to eccentricity of a binary orbit I Peters & Matthews 1963), the initial distribution of eccentricity affects the final expected 
amplitudes. However, we assume that the orbit becomes circular immediately. 
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where M c hi rp = [MiM 2 (Mi + M 2 )~ 1 ^ 3 ] 3 ^ 5 is the chirp mass of the system and c is the speed of light. f p is the reciprocal 
of the proper rest-frame period of the binary, which is related to the observed frequency /, of the gravitational wave 
from the binary in a circular orbit as / = 2/ p /(l + z). D(z) is the comoving distance to the binary 

c f z dz' 

D{z) = tt / . 6 

Ho Jo A /n m (i + z'f + n A 

Next, we calculate v(M\,M 2 ,z) dM\ dM 2 dz, the number of coalescing SMBH binaries with mass Mi and M 2 in 
mass intervals dM\ and dMi per observers' time occurring at z in a given redshift interval dz, 

s dV , 
v{M x ,M 2 ,z) = — n c (Mi ,M 2 ,z) 
ato 

' dV -n c {M x ,M 2 ,z), (7) 



(1 + z) dt p 



where to is the observers' time, t p is the proper rest-frame time at z and dV is the comoving volume element at 
z ~ z + dz given by 

^^^{l + z)cD{zf. (8) 



The timescale emitting gravitational waves of a binary measured in the rest-frame is 

dt 
'df p 



tgw^U^- (9) 



Since the tim escale in the observer-fram e is rcw.obs = tgw(1 + z ), for a binary in a circular orbit at redshift z, it is 
expressed as ijPeters fc Matthewslll963|) 



r G ^ Ml ,M 2 , = ! (^-) [2nf p }- 8/3 (1 
= 1.2 x 10 4 



/ Mchirp \ 5/3 


[ / 1 


V10 8 MoJ 


10" 7 Hz 



-8/3 

(I + Z)- 5 / 3 yr. (10) 

As a SMBH binary evolves with time, the frequency becomes higher. We assume that the binary orbit is quasi- 
stationary (i.e. phase evolution time scale is less than tqw) until the radius equals to 3i?s, where i?s is the Schwarzschild 
radius : the radius of the innermost stable circular orbit (ISCO) for a particle and a non-rotating black hole. Then 
the maximum frequency / max is 

l M 2 y/ 2 



= 4.4xl0- 5 (l + ^)- 1 1 + —) Hz, (11) 

where Mi and M 2 are SMBH masses (Mi > M 2 ). The number of coalescencing SMBH binaries as a source of 
gravitational wave source with observed frequency / at redshift z is estimated as 

v(Mi,M 2 ,z)dz dMi dM 2 T G w,obs 0(/ max - /), (12) 

where 8(x) is the step function. Therefore, from equations (0, @, l|10l) and l|13[i . we finally obtain the spectrum of 
the gravitational wave background radiation, 

hl{f) =Jdz dMi dM 2 h 2 s v(Mi,M 2 , z) T GW , obs 0(f max - /). 

dz dMl dM 2 (^^y\ n f)-^(l + z )-^n c (Mi,M 2 ,z)e(f max -f). 

(13) 

We note that the equation 11131) agrees with the result of iPhinnevI ()2001|) except for the step function, 0(/ max — /)• 
The original formulation of Jaffe & Backer (2003) does not include this effect, either. 

As shown in figure the spectrum changes its slope at / ~ l^iHz owing to lack of power asso ciated with th e uppe r 
limit frequency, /max- This feature is consistent with the results of iWvithe fc Loebl l|2003fl and ISesana et alJ ([2004). 
The predicted strain spectrum is h r (f) ~ 10~ 16 (,f /l/iHz)~ 2 / 3 for / < 1/iHz, just below the current limit from the 
pulsar timing measurements bv lLommenl (|2002(1 . In our model, we assume that SMBHs coalesce simultaneously when 
their host galaxies merge. Therefore, the efficiency of SMBH coalescence is maximum and the predicted amplitude of 
gravitational wave spectrum should be interpreted as the upper limit. 
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In figure |3 (a), we plot the spectra from binaries in different redshift intervals. This figure shows that the total 
spectrum of background radiation comes from coalescing SMBH binaries at low redshift, < z < 1. In figure 01 
we plot the SMBH coalescence rate in observers' time unit a year, v[z) = J v{M\, M2, z) dM\ dMi- This figure 
indicates that the coalescence rate at low redshift is lower than the coalescence rate at high redshift. However, the 
main contribution to the background radiation is gravitational wave from the coalescing SMBH binaries at low redshift, 
< z < 1, because the distance from an observer to the SMBH binaries at low redshift is shorter and the mass of 
SMBHs at low redshift is higher. In figure EJb), we also plot the spectra from binaries in different total mass intervals 
(Mtot = Mi + M%). This figure shows that for / > 10 _4 Hz the total spectrum of background radiation comes from 
coalescing SMBH binaries with total mass M tot < 1O 8 M . 

When SMBHs are spinning and/or when their masses are comparable, the definition of ISCO becomes vague and 
our assumption that the cutoff frequency / max corresponds to 3i?s may not be correct. To see the effects of the 
cutoff frequency, we plot the spectra for different values of / max , corresponding to 3i?s, 30i?s, and no frequency cut 
off, respectively (figure (a)). Lowering / max causes a suppression in the stochastic background at high frequencies 
/ < 10~ 7 Hz since a large portion of high frequency modes are cut off. As shown in this figure, the cut off radius does 
not affect the amplitude at the frequency of the pulsar timing measurement / ~ 10~ 9 Hz. 

Our galaxy formation model incorporates dynamical friction and random collision as galaxy merging mechanisms. In 
order to examine the effect of these two galaxy merger mechanisms on the spectrum of gravitational wave background 
radiation, in figure[S] (b), we also plot the spectrum of background radiation of other two models: no rc model and no 
df model. The no df model can not produce higher mass SMBH (see fig. ED>)). Consequently, the spectrum in no df 
model bends at larger frequency since the number of SMBHs with smaller / max decreases. Furthermore, in the no df 
model, the amplitude of gravitational waves from each binary becomes smaller . However, the coalescence rate in no 
df model is higher than the rate in no rc model as shown in figure As a result, the amplitude of the spectrum in no 
df model is roughly equal to the amplitude in no rc model. 



3.2. Gravitational wave burst from SMBH coalescence 

After an inspiraling phase, SMBHs plunge into a final death inspiral and merge to form a single black hole. We 
call a set of a plunge and a subsequent early non-linear ring-down phase as a burst 1 . In this subsection, we esti- 
mate the expected burst e vent rate per observers' time using the amplitude of burst gravitational wave given by 
iTh orne fe; Braginskvl l)1976|) . and the SMBH coalescence rate calculated by our SA model, n c (Mi, M2, z). 

The gravitational wave amplitude is given by 



ti 2 = 



2GF G 



w 



7TC 3 / 2 



(14) 



where Few is the energy flux of the gravitational wave at the observer and f c is the observed characteristic frequency. 
f c from the gravitational wave burst occurring at z is 



3 3 / 2 GM t ot (1 + z) 
M tot 



= 3.9 x 10" 



10 8 Mr. 



(l + z)- 1 Hz, 



(15) 



where M tot is total mass of the black holes. The energy flux from the burst gravitational wave occurring at z is given 

by 



Fgw = 



eM tot c 2 f c 
4nD(z) 2 (l + z)' 



(16) 



where e is the efficiency of the energy release and D{z) is the distance to the source from the observer given by equation 
Jfjjl. From equations 1|14|) . 1|15|) and (|16fl . the characteristic amplitude of the burst gravitational wave is 



''burst " 



3 3 /V/ 2 GAfto t 

2 1 / 2 7TC 2 D(z) 

= 7.8 x 1(T 16 ( —) 
\0.lJ 



V2/ Aft, 



V 10 8 Mq 



1 Gpc 



(17) 



Equations l|15|) and (|17f) show that we can estimate the amplitude and the characteristic frequency of the burst 
gravitational wave, once we know the total mass and redshift of coalescing SMBHs. From ea. I|15fl . ea. (|17|l and 
n c (Mi, M2, z), we obtain nburst(^burst) fc, z) d/iburst df c dz , which is the comoving number density of gravitational 
wave burst events occurring at z in a given redshift interval dz with amplitude /iburst ~ /iburst + rf^-burst and with 
characteristic frequency f c ~ / c + df c . Then, the expected event rates of gravitational wave bursts per observers' time 



7 Later ring-down phase at linear perturbation regime is not included. 
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with amplitude /iburst ~ h hul - st + dh hul - st and characteristic frequency f c ~ f c + df c , tw-st (/iburst, /c) dh hurst df c is given 

by 

nhmst(hburst, fc, z)——dz. (18) 

The integrated event rates of gravitational wave bursts per observers' time with amplitude /iburst ~ /iburst + c//iburst, 
^WstC^burst) dh hurst , are given respectively as follows, 

^burst (/iburst) — / ^burst 

(h hurst , f c )df c . (19) 

In figureEl we plot the total integrated event rates of gravitational wave bursts and integrated event rates in different 
redshift intervals (figure Eta)) and in different total mass intervals in (figure Efb) ). Here we set the e f ficienc y of the 
energy release e — 0.1, while the precise value of this parameter is unknown. iFlanagan fc Hughes! l|1998fl argued 
that the efficiency could reach 10 %, which depends on align ments of spin and parameter choices. In most typical 
events, a conversion efficiency will probably be a few percent ijBaker et al.ll200l|) . From equation l|17|) . one can see 
that the change of efficiency results in the parallel displacement in the horizontal direction in figure H3 The shape of 
z^burstCibm-st) reflects the black hole mass functions and the SMBH coalescence rates, which depend on the complex 
galaxy formation processes. In figure E] (a), one can notice that there are two peaks in the event rate in terms of 
/iburst- A peak at /iburst ~ 10 -17 corresponds to bursts from SMBH binaries with M tot < 10 6 Af Q whose total number 
is the largest at high redshift z > 3. Another peak at /iburst ~ 10 -15 corresponds to bursts from SMBH binaries with 
10 7 M Q < M tot < 10 8 M© whose coalescence probability is the largest at high redshift z > 3. Figure El (b) indicates 
that burst signals with large amplitude (/iburst ^ 10 -15 ) correspond to coalescence of "massive" SMBH binaries with 
Mtot ^ 10 8 M Q occurring at z < 1. This is because the distance from the earth to SMBHs at low redshift is shorter 
and the mass of SMBHs at low redshift is larger. On the other hand, burst signals with amplitude (/iburst ^ 10 -15 ) 
corresponds to coalescence events of "less massive" SMBH binaries with M to t IS 10 7 M Q occurring at z > 2. These 
events dominate the expected burst event rate provided that the sensitivity of the detector is sufficiently good. This 
feature is quite important because it breaks the degeneracy between mass and distance. Our model predicts that the 
expected rates of the coalescence events owing to SMBH binaries with small mass M to t < 10 6 M© at low redshift and 
those with large mass M tQt > 10 8 M Q at high redshift z > 3 are small. 

Figure [3 shows that the expected region for signal of gravitational wave bursts and the instrumental noise threshold 
for LISA, /linst- For r andomly oriented sources, a s ensitivity for a search of gravitational wave bursts in an observation 
time T ob s is given by <|Thornell 19871 lHaehneltlll994|) 

hLt~(j^) 'lOSb/c, (20) 

wher e Sh is the spectral instrumental noise density. We compute h- lns t from the fitting formula for Sh of LISA ijHughesl 
2002). The expected region for Z4>urst [log (/iburst)) log(/ c )] > 1 yr -1 is above this instrumental noise threshold. For 
comparison, we show the region for i<Wst[log(/iburst),log(/ c )] > 1/5 yr _1 and iburst [log ( /iburst ), log (/ c )] > 3 yr _1 in 
figured (b). 

From figureE]and[71 we conclude that the LISA can detect intense bursts of gravitational waves at a rate of ~ l.Oyr -1 
assuming that dominant part of these burst events occur at z > 2. Even in the case of e = 0.001, the LISA can detect 
intense bursts of gravitational waves in one year observation, since /iburst oc e 1 / 2 . In addition, we find that large 
amplitude /iburst ~ 10~ 15 signals correspond to coalescence events of massive SMBH binaries M to t ~ 1O 8 M at low 
redshift z < 1 and small amplitude /iburst ~ 10 -17 signals correspond to less massive SMBH binaries M tot ~ 10 6 M Q 
at high redshift z > 3. 

I! :-. :; on a SA-model l|Kauffmann fc Haehnehll2000|) . iHaehneltl l)2003t) concluded that LISA might detect SMBH 
coalescence events at a rate 0.1 ~ 1.0 yr _1 over the redshift range of < z < 5 although no explicit calculation on 
gravitational wave emission is done. Our result is consistent with his result. 

4. SUMMARY AND CONCLUSIONS 

In this paper, we have estimated the coalescence rate of supermassive black hole (SMBH) bina ries in the c enters 
of galaxies using a new semi-analytic model of galaxy and quasar formation (SA model) given bv lEnoki et ail C2003) 
based on the hierarchical structure formation scenario. Then, we calculated th e spectrum of the gra vitational wave 
background from inspiraling SMBH binaries based on the formulation given by l.laffe fc Backer! l|2003f) and estimated 
the expected amplitudes and event rates of intense bursts of gravitational waves from coalescing SMBH binaries. 

Our SA model includes dynamical friction and random collision as galaxy merging mechanisms, and assumes that 
a SMBH is fueled by accretion of cold gas during a major merger of galaxies leading to a spheroid formation, and 
that SMBHs coalesce simultaneously when host galaxies merge. Many previous other studies have paid attention 
to only SMBH growth and did not take galaxy formation processes into account. For investigating the relations 
betwe en SMBH growth and galaxy formation processes, SA me thods of galaxy an d SMBH formation are suitable 

ineldRMiri !Catlanedl200lt lEnoki etal]l20?i3 iMenci et al.ll20T)3l lOranato et alJl200l . Our 



tools ijKauffmann fc Haehr 

SA model can reproduce not only observational prop erties of galaxies, but also the present SM BH mass function and 
the quasar luminosity functions at different redshifts ijNagashima etalll200ll TEnoki et al JI20?)1 . 
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We have found that the gravitational wave background radiation spectrum for / < 1/iHz has a characteristic strain 
h c {f) ~ 10 _16 (//l^Hz)~ 2 / 3 just below the detection limit from the current measurements of the pulsar timing. The 
slope of the spectrum for / > 1/iHz gets steep owing to the upper limit in frequency set by the radius of the innermost 
stable circular orbit. The stochastic background radiation mainly comes from inspiraling SMBH binaries at < z < 1. 
Therefore, the background radiation can probe inspiraling SMBH binaries at low redshift. 

We have also found that LISA might detect intense bursts of gravitational waves owing to the SMBH coalescence 
events at a rate 0.1 ~ 1.0 yr _1 and that the main contribution to the event rate comes from SMBH binary coalescence 
at high redshift z > 2. Our model predicts that burst signals with a large amplitude correspond to coalescence of 
large mass SMBH binaries at low redshift while those with a small amplitude correspond to coalescence of small mass 
SMBH binaries at high redshift. This prediction can be tested by f uture measurem ents of the amplitude and the phase 
evolution in gravitational waves from inspiraling SMBH binaries fHughej 1:2002). Comparing these predictions with 
observations in future, we can put a stringent constraint on SMBH formation and evolution models. 
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astronomy and cosmology at Durham. NS is supported by Japanese Grant-in- Aid for Science Research Fund of the 
Ministry of Education, No. 14340290. Numerical computations in this work were partly carried out at the Astronomical 
Data Analysis Center of the National Astronomical Observatory Japan. 

REFERENCES 



Armitage, P. J. & Natarajan, P. 2002, ApJ, 567, L9 
Bahcall, J. N., Kirhakos, S., Saxe, D. H., & Schneider, D. P. 1997, 
ApJ, 479, 642 

Baker, J., Brugmann, B., Campanelli, M., Lousto, O. & Takahashi, 

R. 2001, Phys. Rev. Lett., 87, 121103 
Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 
Begelman, M. C, Blandford, R. D., & Rees, M. J. 1980, Nature, 

287, 307 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics, Princeton, 

NJ, Princeton Univ. Press 
Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 
Bond, J. R., Cole, S. , Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 

440 

Bower, R. J. 1991, MNRAS, 248, 332 
Cattaneo, A. 2001, MNRAS, 324, 128 

Cole, S., Aragon-Salamanca, A., Frenk, C. S., Navarro, J. F., & 

Zepf, S. E. 1994, MNRAS, 271, 781 
Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 

319, 168 

Detweiler, S. 1979, ApJ, 234, 1100 

Dunlop, J. S. et al., 2003, MNRAS, 340, 1950 

Enoki, M., Nagashima, M. & Gouda, N. , 2003, PAS J, 55, 133 

Ferrarese, L. & Merritt, D., 2000, ApJ, 539, L9 

Flanagan, E. E., & Hughes, S. A. 1998, Phys. Rev. D, 57, 4535 

Fukushige, T., Ebisuzaki, T. & Makino, J., 1992, ApJ, 396, L61 

Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., 

Filippenko, A. V., Green, R., Grillmair, C, et al., 2000, ApJ, 

539, L13 

Gould, A. & Rix, H. 2000, ApJ, 532, L29 

Granato, G. L., De Zottie, G., Silva, L., Bressan, A., & Danese, L., 

2004, ApJ, 600, 580 
Haehnelt, M. G., 1994, MNRAS, 269, 199 
Haehnelt, M. G., 2003, Class. Quantum Grav., 20, S31 

Hughes, S. A. 2002, MNRAS, 331, 805 

Islam, R. R., Taylor, J. E. & Silk, J., 2003, astro-ph/0309559 

Jaffe, A. H. & Backer, D. C, 2003, ApJ, 5S37TTr5 

Kauffmann, G., White, S. D. M., & Guiderdoni, 1993, MNRAS, 

264, 201 

Kauffmann, G., & Haehnelt, M. G. 2000, MNRAS, 311, 576 

Kauffmann, G., et al., 2003, MNRAS, 346, 1055 

Kodama, T., & Arimoto, N. 1997, A&A, 320, 41 

Kormendy, J., & Richstone, D., 1995, ARA&A, 33, 581 

Lacey, C. G., k. Cole, S. 1993, MNRAS, 262, 627 

Lommen, A. N. 2002, in the Proceedings of the 270. WE-Heraeus 
Seminar on Neutron Stars, Pulsars and Supernova Remnants, 
2002, Physikzentrum Bad Honnef, eds. W. Becker, H. Lesch & 
J. Truemper. 



Magorrian, J., et al., 1998, AJ, 115, 2285 

Makino, J. 1997, ApJ, 478, 58 

Makino, J., & Hut, P. 1997, ApJ, 481, 83 

McLure, R. J., Kukula, M. J., Dunlop, J. S., Baum, S. A., O'Dea, 

C. P., & Hughes, D. H. 1999, MNRAS, 308, 377 
Menci N., Cavaliere A., Fontana A., Giallongo E., Poli F., & 

Vittorini V., 2003, ApJ, 587, L63 
Menou, K., Haiman, Z., & Narayanan, V. K. 2001, ApJ, 558, 535 
Merritt, D., & Ferrarese, L. 2001, ApJ, 547, 140 
Merritt, D., & Ferrarese, L. 2001, MNRAS, 320, L30 
Mihos, J. C, & Hernquist, L. 1994, ApJ, 431, L9 
Mihos, J. C, & Hernquist, L. 1996, ApJ, 464, 641 
Milosavljevic, M. & Merritt, D., 2003, ApJ, 596, 860 
Nagashima, M., Totani, T., Gouda, N., & Yoshii, Y. 2001, ApJ, 
557, 505 

Nagashima, M., Yoshii, Y., Totani, T., & Gouda, N. 2002, ApJ, 
578, 675 

Peters, P. & Matthews, J. 1963, Phys. Rev., 13, 435 
Phinney, E. S., 2001, astro-ph/0108028 

Portegies Zwart, S. F., Baumgardt, H., Hut, P. & Makino, J. 2004, 

Nature, 428, 724 
Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 
Rajagopal, M. & Romani, R., 1995, ApJ, 446, 543 
Salucci, P., Szuszkiewicz, E., Monaco, P., & Danese, L. 1999, 

MNRAS, 307, 637 
S esana, A., Haardt , F., Madau, P., & Volonteri, M. 2004, 

astro-ph/0401543 
Somerville, K. S., &', Kolatt, T. S. 1999. MNRAS. 305. 1 
Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 
Thoul, A. A., & Weinberg, D. H. 1996, ApJ, 465, 608 
Thorne, K. S. & Braginsky, V. B., 1976, ApJ, 204, LI 
Thorne, K. S. 1987, in Three Hundred Years of Gravitation, ed. 

S. W. Hawking & W. Israel (Cambridge: Cambridge University 

Press) 

Tremaine, S. et al., 2002, ApJ, 574 ,740 

Volonteri, M., Haardt, F. & Madau, P., 2003 ApJ, 582,559 
Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 590, 691 
Yu, Q. 2002, MNRAS, 331, 935 



Gravitational Wave from Supermassive Black Hole Coalescence 



TABLE 1 

Model parameters. 



Cosmological parameters 


Astrophysical parameters 


no n -7 n -7 n n n nnL-!l 


V hot (km s _i ) o hot t% (Gyr) a, / bu i KC T 




Fig. 1. — Averaged SMBH mass growth rate, (Mbh), of the model. The solid, dot-dashed and short dashed lines indicate SMBH mass 
growth rate of total, due to SMBH coalescence and due to gas accretion, respectively. 
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log[M BH /M SUN ] log[M BH /M SUN ] 



Fig. 2. — (a) Black hole mass function of the model at a series of redshifts. The solid, dot-dot-dashed, dashed and dot-dashed lines 
indica te the results at z = 0, 1,2, and 3, respectively. The symbols with errorbars are the present black hole mass function obtained by 
ISalucci et all tl99flh . (b) Black hole mass function of models for three models at z = 0. The solid, dot-dashed and short dashed lines 
indicate both model, no random collision model (no rc model) and no dynamical friction model (no df model), respectively. Both model 
includes dynamical friction and random collision as the galaxy merging mechanism. 
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log(f [Hz]) log(f [Hz]) 

Fig. 3. — Spectrum of gravitational wave background radiation, h c (f), from SMBH binaries in different redshift intervals (a) and in 
different total mass ranges(b). (a) The total spectrum (solid line) and the other lines show those in different redshift intervals < z < f 
(dashed line), f < z < 2 (dot-dashed), 2 < z < 3 (shor t dashed) and 3 < z (long dashed). The filled reverse triangle shows the current 
limit from pulsar timing measurements ILom mcn (2003). 3 ) The total spectrum (solid line) and the other lines show those in different 
total mass intervals M to t < f0 6 M Q (dashed line), 1O 6 M < Af to t < 10 7 M Q (dot-dot-dashed), f0 7 M Q < M to t < 10 8 M Q (long dashed), 
f0 8 M Q < Mtot < 1O 9 M (dot-dashed) and 1O 9 M < M tot (short dashed). 




z 

Fig. 4. — SMBH coalescence rate in observers' time unit a year, v(z). The solid, dot-dashed and short dashed lines indicate both model, 
no random collision model (no rc model) and no dynamical friction model (no df model), respectively. Both model includes dynamical 
friction and random collision as the galaxy merging mechanism. 
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log(f [Hz]) log(f[Hz]) 

FlG. 5. — Spectrum of gravitational wave background radiation, h c (f), in differen t / ma x (a ) and i n different galaxy merger models (b). 
The filled reverse triangle shows the current limit from pulsar timing measurements Lommcn (2002). (a) The solid, dot-dashed and short 
dashed lines indicate the results with /max corresponding to 3i?s, 30i?s an d no frequency cut off, respectively, (b) The solid, dot-dashed 
and short dashed lines indicate both, no random collision model (no rc model) and no dynamical friction model (no df model), respectively. 
Both model includes dynamical friction and random collision as the galaxy merging mechanism. 
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log (h burst ) log (h burst ) 

Fig. 6. — Integrated event rate of gravitational wave burst per observers' time unit a year fburst [log(/iburst)]- ( a ) The total integrated 
event rate (solid line) and the other lines show those in different redshift intervals < z < 1 (dot), 1 < z < 2 (dot-dashed), 2 < z < 3 (short 
dashed) and 3 < z (dot-dot-dashed), (b) The total integrated event rate (solid line) and the other lines show those in different black hole 
mass intervals M to t < 10 6 M Q (dashed line), 10 6 M Q < M to t < 10 7 Mq (dot-dashed), 10 7 M Q < M to t < f0 8 M Q (dot) and I0 8 M Q < M tot 
(dot-dot-dashed) . 




log(f c [Hz]) log(fc [Hz]) 

Fig. 7. — Expected signals of gravitational burst from SMBH. (a) The horizontally hatched area shows the region, 
^burst [log(/i b urst)T°g(/c)] > ly r 1 - The solid curve indicates the instrumental noise threshold for one year of LISA observations, (b) 
The vertically hatched area shows the region, fburst [l°g(^burst)) l°g(/c)] > l/5yr _1 and the diagnal cross-hatched area show the region, 
i/b urs t[log(/ib urat ), log(/ c )] > 3yr _1 . The dot-dashed and the short dashed lines indicates the instrumental noise threshold for 5 year and 
1 /3 year of LISA observations, respectively. 



